NPS55-85-031 



NAVAL POSTGRADUATE SCHOOL 



Monterey, California 



TRANSFORMATION OF DISPLACEMENT 
AND ORIENTATION CORRECTION DATA 

GENE GYGAX 
ROBERT R. READ 

DECEMBER 1985 



Approved for public release; distribution unlimited. 
Prepared for: 

Naval Undersea Warfare Engineering Station 
3 voort, WA 98345 




TECHNICAL 



FedDocs 



D 208.14/2 
NPS-55-85-031 




NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 



Rear Admiral R. H. Shumaker 
Superintendent 



D. A. Schrady 
Provost 



Reproduction of all or part of this report is authorized. 



This report was prepared for Naval Undersea Warfare Engineering 
Station by: 



Unci ass f fied 



igcuftirv CLASS. ^CATION QP This PAGE 



Of 1 



REPORT DOCUMENTATION PAGE 



a. Rt-'ORT SECURITY CLASSIFICATION 

Unclassi fied 



lb RESTRICTIVE MARKINGS 



a. SECURITY CLASSIFICATION AUTHORITY 



D. DECLASSIFICATION /DOWNGRADING SCHEDULE 



3 distribution/ availability of report 

Approved for public release; distribution 
is unlimited 



PERFORMING ORGANIZATION REPORT NUMBER(S) 

NPS55-85-031 



S MONITORING ORGANIZATION REPORT NUMBERS) 



a. NAME OF PERFORMING ORGANIZATION 

Naval Postgraduate School 



6b OFFICE SYMBOL 
(If applicable) 



7 a NAME OF MONITORING ORGANIZATION 

Naval Undersea Warfare Engineering Station 



c. ADDRESS (Cry, Srate, and ZiPCode) 

Monterey, California 93943-5100 



7b ADDRESS (C/ry, Stare, and ZiPCode) 

Keyport, Washington 98345 



a. NAME OF FUNDING /SPONSORING 
ORGANIZATION 

aval Undersea Warfare Engineerii< 



c.^auWIsS (City, Srare. and ZIP Code) 



So OFFICE SYMBOL 
(If applicable) 

*3 



9 PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER 

'IIP O/H 



10 SOURCE Of FUNDING NUMBERS 



Wa 98345 


PROGRAM 


PROJECT 


TASK 


//OR*. ^NiT 




ELEMENT NO 


NO 


NO 


ACCESSION NO 

N00253 













1 TITLE {Include Security Classification) 

TRANSFORMATION OF DISPLACEMENT AND ORIENTATION CORRECTION DATA 



2 PERSONAL AUTHOR(S) 

Gyqax, Gene and Read, 


Robert R. 






3a. TYPE OF REPORT 

Techni cal 


13d TIME covered 
FROM TO 


14 DATE OF REPORT (Year, Month. Day) 

85, December 


15 PAGE COUNT 


6 SUPPLEMENTARY NOTATION 






- 



7 COSATI CODES 


FIELD 


GROUP | SUB-GROUP 








“r 



18 SUBJECT TERMS {Continue on reverse if necessary ana identity oy diock numoer) 

Bivariate Data, Transformations, Principal Components 



9 ABSTRACT (Continue on reverse if necessary and identify by block numoer) 

The simulation of an underwater vehicle track has been used to develop reference 
distributions for displacement and orientation corrections to the test range sensor 
arrays. These distributions show unusual behavior. Of particular concern are the 
scatter plots of the magnitude of displacement and orientation rotation angle, which 
exhibit a "hairpin" shape. The mechanisms producing this shape are interpreted and 
data trans formations are developed that produce scatter plots of more regular form. 
In addition, some plotting methods are applied to provide some visual information 
about the joint distribution of the transformed data. 



0 DISTRIBUTION/ AVAILABILITY OF ABSTRACT 
E UNCLASSIFIED/UNLIMITED □ SAME AS RPT □ DT'C USERS 



21 ABSTRACT SECURITY CLASSIFICATION 



2a NAME OF RESPONSIBLE INDIVIDUAL 

Robert R. Read 



2 2d TELEPHONE (include AreaCoae) 

(408)646-2382 



22 c. OFFICE SYMBOL 

Code 55Re 



r~\ r r> r* • • « _ ( 



APR pQition ^<iv ne used until pxnausted 

A I t a * O rs jro jO'.P'^Te 



CLASS. “*CA T 



T- c D AGE 



I. Introduction 



Imagine the scatter plot of data from two strongly correlated variables 
with the provision that only the magnitudes of the individual variables are 
recorded. Such data can produce "hairpin" shaped scatter plots (see Figure 
1) and it would be useful to develop methodology for straightening the 
hairpin, i.e. for assigning negative signs to some of the values. Then 
perhaps one or more of the following could be achieved with the transformed 
data: 

i. The scatter plots could be fitted with bivariate density 
functions, whose contours are convex sets, 
ii. The scatter plots could be fitted with bivariate normal 
density functions. 

In general this problem is a difficult one. Some, but not all, of the 
points that fall near a coordinate axis need to be reflected across that 
axis. A rule for making such selections should have some depth. The 
present work does not treat the general problem. We are more interested in 
expediency so that the general nature of the straightened scatter plot may 
be better understood. Also, our data structure contains concomitant 
information that will be exploited. 

The mother project, of which this is a part, deals with the maintenance 
of calibrations of the 3~D underwater test ranges operated by NUWES. The 
tracking of vehicles is performed by a grid work of three dimensional sonar 
arrays located underwater at the bottom. A number of these individual 
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arrays contribute tracking information which is integrated to reconstruct 
the path of the targeted vehicle. Mismatches occur when the splicing of 
individual contributions is done and the possibility that they are caused by 
sensor movement must be faced. 

Fuller details can be found in References 1 , 2, 3. Of immediate 
concern is the product of the program KEYMAIN which takes data from the 
array overlap region and estimates the relative displacement and orientation 
corrections that minimize the track splicing mismatch. To judge the 
significance of these estimates one needs a reference distribution which 
describes probabilistically the amount of variability that must be tolerated 
in a given situation. A simulation model, whose purpose is to generate such 

reference distributions, has been developed and is described in Reference 1. 

It simulates the path of the vehicle, complete with inherent or random 

variability, for each of the two sensor arrays in an overlap region. The 

results are passed to the program KEYMAIN, which estimates the corrections 
to one array relative to the other one. 

The output of KEYMAIN contains a pair of three dimensional vectors. 

The first is the set of three displacement corrections (downrange, 
crossrange, and vertical) for an array T s location. The second is the set of 
three Euler angles that describe the orientation correction for the array. 
Use is made of two other ways to describe these latter corrections: 

(i) The orthonormal matrix B which provides the information in the form of 
a linear transformation. (ii) The axis direction numbers and the angle of 
rotation. 



2 



It is hoped that these six pieces of information can be replaced with 
two, namely the magnitude of the displacement and the angle of rotation. 

The joint distribution of these two can be economically plotted on a piece 
of paper, and it seems reasonable that this should be adequate to support 
the decision processes mentioned above. Such are the data in Figures 1 and 
2 . 



Some additional unusual behavior is exhibited in Figure 2; namely, 
there seems to be a "base line of support" which the bivariate data do not 
cross. The mechanism that produces this is not understood at this writing, 
and this feature will not be treated. 

The organization of the report is as follows: Section II contains the 

mathematical development of methodology to " straighten the hairpins" and 
perform ad hoc graphic tests for the bivariate normal distribution. Section 
III contains the application of this methodology to our six cases of 
simulated data. A summary follows. 

II . Mathematical Support 

The simulated data are six dimensional and are best viewed as two 
strongly correlated 3"D sets; the three components of displacement and the 
three components of orientation. As a data reduction technique let us draw 
down to two dimensions; the magnitude of displacement d, and- the magnitude 
of the angle, 0, of reorientation about an axis of rotation. It is the 
scatter of 0 on d that exhibits the unusual behavior (e.g. Figure 1). 
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The hairpin shape is readily interpreted once one is aware that 
displacement and orientation are strongly correlated, and that the 3~D data 
swarm of displacement is very long and narrow. That is, the variability of 
displacement is dominated by its first principal component. At one end of 
this principal component both 0 and d are large. They both become smaller 
as one moves away from the extreme and toward the center of the principal 
component. But one of these, say 0, reaches its minimum first and starts 
becoming larger again before the other, d, achieves its minimum value. This 
may cause an observable "arch on the pin". Then the data fills out the 
other leg of the pin as we move to the other extreme of the principal 
component. Figures 3 and bear this out. These are scatter plots of d, 
resp. 0, on the first principal component, u. The key to straightening the 

hairpins may be found in assigning negative values to the left leg of these 

"v" shaped scatter plots. 

This section contains the mathematical support needed to straighten the 
hairpin and perform an ad hoc test for bivariate normal distribution to the 
result. It contains three subsections: (i) Development of formulas 

relating the three ways to represent the orientation corrections. (ii) 
Discussion of methods to assign negative signs to some of the d and 0. 

(iii) Development of an ad hoc test for bivariate normal distribution. 

They are presented in order. 






Subsection (i). Three Ways to Represent the Orientation Correction 



The three ways are: 

1) Orthonormal matrix (B) . 

2) Euler angles (<(>!, <J> 2> <t> 3 ). 

3) Axis (x) and rotation (0). 

Each has its advantages. We proceed to describe them and show how to 
transform from one to the other. 

T 

1. An orthonormal matrix B is defined by the property B B = I (the identity 
matrix) and may be viewed as a linear transf ormat ion that preserves 
distance. Thus, if x and y are two vectors, the distance separating them 
can be computed from 

| |x - y| | 2 = (x - y ) T (x - y) - E(x. - y^ 2 . 

To illustrate the distance preserving property, let x f = Bx and y* = By 
represent the images of x and y under the linear transf ormation B. Then 

| | x 1 - y * I I 2 = I I Bx - By | | 2 = | | B ( x - y ) | | 2 = (B (x - y) T (B(x - y)) 

= (x - y) T B T B(x - y) = (x - y) T (x - y) = ||x - y | [ 2 

Thus the distance between x and y and the distance between their images is 

the same. These tranformat ions are the rigid rotations. Because of the 
T 

constraint, B B = I, the matrix B has only three degrees of freedom. Three 
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are lost because each column must have length unity, and three more are lost 
because each pair of columns must be orthogonal. 

2. The three Euler angles specify the rotations that take place in the base 
planes of the coordinate system as each axis is held fixed. Since the 
positions of these planes change with each rotation, it is necessary to 
specify the successive order in which the axes are held fixed. We choose to 
first hold the x-axis (downrange) fixed and rotate the y - z plane through 
all angles of <$> x ; then hold the y-axis (crossrange) fixed in its new 
position so as to rotate the (new) x - z plane through an angle <j> 2 ; followed 
by holding the z-axis (vertical) fixed in its (twice) tranformed position 
and rotate the x - y plane through an angle <f> 3 . 

In navigation work, these three angles are called, roll, pitch, and 
yaw. In the NUWES software they are designated as Y-TILT, X-TILT, and 
Z-ROT. 

The Euler angles are easily related to the B-Matrix. Each is a 
rotation in a plane of two coordinate axes while the third is held fixed. 

It follows that each rotation has its own repr esentation as an orthonormal 
matrix. Thus, letting 

S = sin( 4> ) and C = cos (((>.) 
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for i = 1, 2, 3, we can specify: 




(Note that we have specified clockwise, or "back" rotation. Then, using the 
stated order of application, 



B = 83828! 



and the algebraic result of these multiplications is 





( C 2 C 3 


^ \ S 2 & 3 + 


-(C l S 2 C 3 - SjSj) 


B = 


j C 2 S 3 


-(S 1 S 2 S 3 - C 1 C 3 ) 


-(C,S 2 S 3 + S,C 3 ) 




ls 2 


S 1 C 2 


CjC 2 



Given B, the Euler angles can be recovered in the following way. 
Compute principle value solution using 



sin((j> 2 ) = b 31 sin($ x ) = b 32 /cos(4> 2 ) sin(<t> 3 ) = b 2 1 /cos(4) 2 ) 



where 



~tt/2 < . < tt/ 2 for i = 1 ,2,3. 

1 
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Then check for possible modification of 4> x and 4> 3 using the signs of and 
b 33 . That is, make replacements, 

<- tt - 4> ! if b 33 <0 
<f>3 <- tt - 4> 3 if b x ! < 0 

3. The axis and rotation representation gives recognition to the fact that 
a rigid turn in three space must involve the turning of the entire space 
through an angle of 0 about an axis x. The only points unmoved by this 
transformation are the points on the axis. Thus 

Bx = x 

and the axis x is identified as the eigenvector of B having eigenvalue one. 
In order to make it unique, let us agree that x is normalized so that 
||x|| 2 = 1 and it lies in the half space defined by a non-negative first 
component, x x ^ 0. Note there are two degrees of freedom in this x. 

To determine the rotational angle 0, let us set up a righthanded 

coordinate system with x as its first axis, and then determine the angle of 

rotation resulting from applying B to any vector orthogonal to x. Our goal 

is to construct a change of basis matrix C which has x as its first row. It 

remains to specify the second and third rows, y and z. Since x is pointed 

T 

in the half space Xj > 0, let us apply the Gram-Schmidt process to e 2 = 
(0,1,0) in order to specify y. 
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I.e. 



y = e 2 - (x j ,e 2 )x 



followed by a normalization to give it length one. Thus 



t 

y = (-XiXa, 1 - X 2 2 > ““X 2 X 3 ) / / l - x 2 2 



In the construction of z we want to be assured that the result is a 
right handed system, so we use the vector cross product. 



z = x x y = det 



e i e 2 e 3 

x ! x 2 x 3 

y 1 y 2 y 3 



where the representation by the determinant is a schematic one; inferring an 
expansion by the first row, whose components are the unit coordinate 
vectors . 



The result is 



z T = (-x 3 , 0, x j ) / 1 - x 2 2 



Let us apply our change of basis matrix C 

T 

C = (x,y,z) 
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to the vector By. From geometrical considerations, we have By rotates y 
through the angle 0 and its representation in the basis C must be 



The third component can be inverted to compute the angle 9. A check of the 
sign of the second component will tell us the quadrant to which 9 belongs. 
Thus x and 9 can be constructed from the orthonormal matrix B. 

To reverse the process, i.e. to construct B from the axis x and angle 9 
we can proceed as follows: 

First construct the orthonormal matrix C from x exactly as before. 

Then apply the matrix 



to C. This will rotate the plane of the second and third coordinates 
through the angle 9 with the first coordinate held fixed. One can verify 
easily that 



CBy = 





B*CBx = e x B* CBy = e 2 



B* CBz = e 3 , 



1 0 



i .e. 



b*cbc t = I 



and from this equation we may solve for B. 

T T 

B = C B* C 



Subsection (ii). Assignment of Signs. 

One first needs the values of u, the first principal component of the 
3~D displacement data swarm. We simply state that if x is a set of 3“D 
displacement vectors, and p is the eigenvector corresponding to the largest 
eigenvalue of the covariance matrice of { x } , then the values of u are the 
dot products 



u = p • (x - x) 

where x is the centroid (average) of the data swarm. 

When one regresses either d or 0 on u, the resulting scatter plot has 
the appearance of a "v", or an absolute value function. A difficulty is 
encountered in that the minimum of d or 0 is not necessarily (sufficiently 
close to) zero. This creates the need to estimate values d^ and 0^ so that 



d-d. 



0 0 



1 1 



the replacements 



will produce a continuous regression scatter plot of d (resp. 0) on u after 
the ordinate values on the left leg have by has been assigned negative 
signs. Our approach to this problem has been an expedient one: First re- 

index ©,d, and u according to the rank of u (i.e. increasing order). Next 
select the minimum d(0) of d(0) and let i(j) be the index of d(0) . Then let 

d Q = median (d._ 1 , d. , d. +1 ) 

o j-r j j+i 

It will be seen that this choice is modestly successful. 

Should the overall tr ansf ormation technique be sufficiently promising, 
it will be necessary to develop a more refined method to estimate d^ , 0^. 

An obvious first thought in this direction would be to fit regression lines 
to the two separate legs and use the intersection point. This seems 
attractive except in those cases for which either d^ or 0^ is naturally 
close to zero. There some of the values in the scatter plot have already 
been reflected, e.g. observe the "mushy" bottom in Figure 3. This situation 
interferes with our ability to fit good lines. 



12 



Subsection (iii) Graphical Method to Judge the Appropriateness of the 
Bivariate Normal Model. 



A two dimensional scatter plot that appears consistent with a straight 
line regression model might be studied further for distributional properties 
to determine whether the bivariate normal model is tenable. An important 
and easily exploitable property of this model is that it can be tranformed 
to a circular normal random vector (X,Y) and, upon transforming to polar 
coordinates 



R 2 = X 2 + Y 2 
PV * tan' 1 (Y/X) 
0 = PV 
= PV + IT 
= PV - IT 



for X > 0 

for X < 0 and Y > 0 
for X < 0 and Y > 0 



we have the exponential (1/2o 2 ) distribution to describe R 2 and the uniform 
distribution to describe 0. The two quantities are independent and 
the two distributions are readily verified with the use of simple plotting 
methods . 

It seems wise to take a moment and verify these properties. The 
circular normal density is 

1 { 1 

f(x,y) « - 2 ^r exp - -^t-(x 2 - y 2 ) 



1 3 



The inverse polar transf ormation is 



x = r cos(O) y = r sin(0) dxdy = rdrd0. 



Use of this substitution produces 



6(r,0)) - re' rV2 ” 2 - 8„(r>^9> 



, v -r 2 /2a 2 n . 

where g D (r) = — T" e , 0 < r < 00 

K 0 

i 

and %ri-©) * ' 2 ~~ > — it < 0 S it. 



The result will be completed when we make the change 



s = r 2 ds = 2rdr 



so that 



/ \ 1 -s/2a 2 ^ 

g s (s) = jp- e > 0 < s < “. 



Thus (H) is distributed uniformly as (— ir,ir) 



and S = R 2 is exponential (1/2o 2 ). 
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Histograms of © and £n(R 2 ) should exhibit a constant value behavior, 
and a scatter plot of © on R(or R 2 ) should not suggest dependence. 

The preliminary step of transforming the original data, say (W,Z) into 
potential circular normal form needs to be described: 

First make the replacements 

w «- (w - w)/S 

w 

z «- (z - z)/S 

z 

where w and z refer to the original center of gravity of the data -and S and 

w 

are the sample standard deviations. Then compute the correlation 
coefficient 



1 N 



where N is the number of data vectors. This done, the joint distribution of 
W, Z is viewed as being a standard bivariate normal with correlation p. The 
final step is to apply the transf ormation 

x = w 

y « (z -pw)// 1 - p 2 . 
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It is easily verified that (X,Y) has a circular normal distribution and the 
variance is unity. 



The additional step, taking (X,Y) to (R 2 ,® ) was described earlier. In 
addition to (or in place of) the histogram plots mentioned earlier, one can 
apply quantile plots, i.e. theoretical quantiles plotted against empirical 
quantiles. The value of these plots lies in the fact that they should 
produce reasonable looking straight lines if the theoretical distributions 
are well chosen. There is great sensitivity in the tails (extremes) of 
these plots however. 

The empirical quantiles are obtained from the order, statistics: 




W = 3 



= i/ (N + 1 ) 



for i 



, • • • 



N 



The corr esponding theoretical quantiles are obtained from 



Sej p i } = ~ 1} 

Qg(P i ) = “4n(1 - p.) 
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III. Analysis of Application. 



The data sets used are those available from Reference 1; a Mark 10 
track at Nanoose on 9 September 1982. The particulars are: 

Point No. of 

Data Set Sensor Pair Count Interval Usable points 



ONE 


4,5 


3933,3999 


67 


TWO 


4,5 


4028, 1)097 


70 


THREE 


4,5 


4405,4458 


53 


FOUR 


5,6 


6909,7001 


93 


FIVE 


5,6 


7079,7162 


82 


SIX 


5,6 


7647,7738 


87 



Each of these six situations were simulated 200 times, and the 
displacement and orientation corrections were computed each time. Scatter 
plots of 0 on d appear in Figure 1 and 2. To get a reasonable view of these 
data, it was necessary to remove some wild simulated points. Specif ically : 



Data Set 


ONE 


TWO 


THREE 


FOUR 


FIVE 


SIX 


Outliers Removed 


1 


24 


41 


0 


0 


0 



The outlier removal rule was: Drop the values if either d > 100 ft or 

0 > 0.5 radian. 
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The hairpin shape is obvious for data sets one and two. It will be 
seen that it is also clearly present in data sets three and six as well. 

Data sets four and five seem to have a boundary. No interpretation for this 
is offered at this time. The hairpin effect may be explained as the scatter 
plot of the absolute values of two variables that are (otherwise) linearly 
related. 

Some other features of Figure 1 and 2 are noted in passing. Generally 
the first three data sets have larger values for d and 0 than do the other 
three. This fact could be related to one or both of the following facts: 

(i) The valid point count size is smaller for the first 
three; 67, 70, 53 vis a vis 93, 82, 87. 

(ii) The sensor array pair is different; (4,5) vis a vis 
(5,6). 

Figures 3 and 4 contain the scatter plots of d and 0 versus the first 
principal component of their respective displacement data sets. The 
efficacy of using a principal component to represent a multivariate data set 
is dependent upon the percent of total variability that is represented by 
that component. These figures are recorded below: 



Data Set 


ONE 


TWO 


THREE 


FOUR 


FIVE 


SIX 


Percent of Total Variance 


97.0 


74.6 


59.3 


85.5 


99.4 


99. 1 


Total Variance 


6349.6 


14564.2 


275935.4 


13.1 


1 19.8 


793.8 
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It is interesting to note that the two ’’non hairpin” cases, four and five, 
are of conspicuously low variability. The appearance of the ”v” shape is 
suggestive of the absolute value of an otherwise signed variable. Since 
data sets 1, 2, 3, and 6 show this shape, they are the ones that can produce 
hairpins in the scatter plots of 0 and d. Two of these four hairpins ( 3 
and 6) have such narrow arches that the scatter plots of the two legs do not 
separate (see Figures 1 and 2). We note that the plot of d vs u for Data 
Sets four and five is suggestive of a hairpin whose bend point is not in the 
center but deep in the tail. The corresponding plots of 0 vs u do not 
confirm the effect however. 

The assignment of negative signs to the left legs of these ’V’ shaped 
plots must be preceded by the translations. 





d <- 


d - d Q 


and 


0 <- 


°0 




Using the local 


median of 


three 


technique 


explained 


in Section II, 


values are 














Data Set 


ONE 


TWO 


THREE 


FOUR 


FIVE 


SIX 


D o 


2.20 


6.07 


15.10 


15.19 


9.74 


15.19 


e o 


.0028 


.0030 


.0023 


.0077 


.001 6 


.0017 



Having made these translations, negative signs are attached to all values of 
d(0) whose indices (computed from the ranks of u) are less than or equal to 
the index of d,(0) = zero (i.e. formerly d o ,(0 Q )). 
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This done, the Figures 1 and 2 are converted to their signed versions, 
Figures 5 and 6. The result appears to be reasonable in the four hairpin 
designated cases (Data Sets 1, 2, 3 and 6). A look at the center of the 
plots for Data Sets 2 and 3 suggests that the values for d Q and 0 Q might be 
better chosen. The result for Data Set 5 may not be as good as the original 
untransf ormed plot. 

Turning to the application of our ad hoc test for the bivariate normal 
distribution , recall that the data are transformed first to the circular 
normal and then to polar form. Figure 7 contains the scatter plots of the 
squared radius vs the polar angle. It was hoped that these could show 
independence. The zero regression aspect of independence seems to be 
present, but at least three of the plots (Data Sets 1, 3, 5) exhibit some 
sort of cyclical effect; large radii are associated with regularly spaced 
polar angles. 

Figures 8 and 9 contain Q - Q plots; empirical quantiles in the 
horizontal and theoretical quantiles in the vertical. To avoid crowding, 
every third point was plotted. In order to support the distributional 
assumptions, all of these plots should be reasonably straight lines. It is 
our opinion that the departures from straight lines are sufficient to cast 
doubt on any precise use of the bivariate normal distribution for the 
scatter plots in Figures 5 and 6. 
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FIGURE 1: ROTATION VS DISPLACEMENT USING SENSORS FOUR AND FIVE 

DATA SET ONE, N= 199 OATA SET TWO, N=176 
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MAGNITUDE OF DISPLACEMENT (FEET) 



FIGURE 2: ROTATION VS DISPLACEMENT USING SENSORS FIVE AND SIX 

DATA SET FOUR, N = 200 DATA SET FIVE. N=200 
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MAGNITUDE OF DISPLACEMENT (FEET) 



FIGURE 3: DISPLACEMENT/ROTATION VS PRINCIPLE COMPONENT 
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ORDERED PRINCIPLE COMPONElfT ORDERED PRINCIPLE COMPONENT ORDERED PRINCIPLE COMPONENT 



FIGURE 4: DISPLACEMENT/ROTATION VS PRINCIPLE COMPONENT 
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ORDERED PRINCIPLE COMPONENT ORDERED PRINCIPLE COMPONENT ORDERED PRINCIPLE COMPONENT 



FIGURE 5: SIGNED ROTATION VS DISPLACEMENT, SENSORS FOUR AND FIVE 

DATA SET ONE, N=199 DATA SET TWO, N=179 
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FIGURE 6: SIGNED ROTATION VS DISPLACEMENT, SENSORS FIVE AND SIX 

DATA SET FOUR. N=200 DATA SET FIVE, N-200 
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FIGURE 7: DATA SCALED AND TRANSFORMED TO POLAR COORDINATES 
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FIGURE 8: QUANTILE PLOTS, SENSORS 4 AND 5 
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IV. Summary 



The methodology developed here is certainly useful but there are 
aspects of the data that are not yet understood. More and higher quality 
data is needed in order to continue the development. The main goal is to 
substitute two dimensional data for six dimensional data. In doing so, it 
is hoped that one can profitably use popular bivariate probability laws to 
describe these references distributions. 
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